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VLASOV EQUATION EIGENVALUES AND EIGENVECTORS FOR 
FOURIER- HERMITE DISPERSION MATRICES OF 
ORDER GREATER THAN 1000 

By Frederick C. Grant 
Langley Research Center 

SUMMARY 

The connection between the Van Kampen and Landau representations of the Vlasov 
equations has been extended to Fourier-Hermite expansions containing more than 
1000 terms by taking advantage of the properties of tridiagonal matrices. These numer- 
ical results are regarded as conclusive indications of the nonuniformly convergent behav- 
ior of the approximation curve in the limit of an infinite number of terms and represent 
an extension of work begun by Grant (1967) and by Grant and Feix (1967). 

INTRODUCTION 

In a fully ionized plasma such as the solar corona, which has a large number of 
electrons in a Debye cube (the edge of which measures the screening distance about an 
introduced charge), the effect of collisions is negligible and the Vlasov equations describe 
the space-time evolution of the single-particle distribution function. 

The nonlinear Vlasov equations have been studied only for special cases, a common 
one of which is assumed herein: a single space and velocity variable, periodic spatial 
boundary conditions, a uniform steady ion background against which the electrons move, 
and an initial Maxwellian electronic distribution function. 

Because of their formidable nature, the Vlasov equations have been treated mostly 
by numerical techniques. However, when ordinary x-v-t space is used, an insuperable 
difficulty arises in the velocity space. Increasingly steep oscillatory gradients with ever- 
shorter wavelengths develop in the distribution function. To circumvent the difficulty, 
transform methods have been introduced, and among these is the Fourier expansion on 
space and the Hermite expansion on velocity which lead, in the nonlinear case, to a set 
of coupled ordinary differential equations for the elements of a matrix, the Fourier- 
Hermite expansion coefficients of the Fourier-Hermite representation. (Each trans- 
formed set of equations is called a representation of the original x-v-t equations.) 



When the nonlinear Fourier-Hermite equations are linearized and an eigenvalue- 
eigenvector solution is sought, another matrix appears, the Fourier-Hermite dispersion 
matrix. At this point contact is made with two classical linearized representations: that 
of Landau (ref. 1) and that of Van Kampen (ref. 2). Landau uses a Fourier- Laplace trans- 
form on space and time, while Van Kampen performs essentially a Fourier transform on 
both space and time. As it happens, the continuum of eigenvalues in Van Kampen' s repre- 
sentation are all undamped (real) while the discrete poles of Landau are all damped (com- 
plex). In particular, a pair of highly excited poles, called Landau poles, dominates the 
motion for long times when the initial conditions are stable. The Fourier-Hermite repre- 
sentation has discrete, real eigenvalues, and hence it seems clearly akin to Van Kampen's 
representation. 

When numerical integration of the Fourier-Hermite equations is performed, a 
numerical instability appears at the highest order Hermite index and spreads toward 
the lower. If a Fokker-Planck (weak- collision or distant-encounter) term is introduced 
into the Fourier-Hermite representation, it damps most strongly the coefficients of high- 
est Hermite order. Macroscopic physical entities such as the electric field depend on 
coefficients of the lowest orders, and thus the selective character of the Fokker-Planck 
term makes it possible, with many Fourier-Hermite coefficients, to quench the instability 
with few collisions and little degradation of the results. At this point the fact that the 
Fourier-Hermite representation bridges the gap between Landau and Van Kampen repre- 
sentations was discovered (ref. 3). 

In reference 3 the transition between Van Kampen and Landau representations of 
the Vlasov equations is shown in terms of a delicate limiting process that requires an 
indefinite increase in the number of terms kept in the approximate Fourier-Hermite dis- 
tribution function and, simultaneously, an indefinite decrease in the strength of the intro- 
duced Fokker-Planck collision term. A full account of the Fourier-Hermite representa- 
tion may be found in reference 4. The full matrix techniques used in reference 3 allowed 
only about 60 terms to be kept because of computer storage limitations. When recursion 
techniques are used in the tridiagonal dispersion matrices, the storage required is pro- 
portional to the matrix order rather than to the order squared. The economy of recur- 
sive methods has allowed the calculation of the eigenvalues and eigenvectors for disper- 
sion matrices larger than 1000 x 1000. 

The ability to treat Fourier-Hermite matrices of orders as great as 1000 carries 
the number of Fourier-Hermite coefficients far enough to make available the appropriate 
value of Fokker-Planck damping for best plasma simulation over the entire practical 
range of expansion coefficient order (several hundred) for the nonlinear problem. 

In the general, nonlinear problem, as in the linear problem, addition of a small 
damping term of proper strength extends the time of good simulation at the cost of slightly 
degraded results. 
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SYMBOLS 


a Q£/3 

B 


C(f) 


E 


e 


F 

F » = dF 
dv 

*o 

f 

f-o 

ft 

i = ^T 
k 

*{} 

m 

N 


Fourier-Hermite expansion coefficient 

strength parameter of Fokker- Planck collision term 

collision term 

determinant of Fourier-Hermite dispersion matrix 
electric field 

magnitude of electronic charge 

exp(-v2/2) 

Maxwellian velocity distribution, 7= — 

\}2v 


Fourier-transform operator 

distribution function of electrons 

Fourier transform of initial perturbation in f 

Fokker-Planck collision term 

Hermite polynomial of order j3 

n 

wave number, — ^ : Boltzmann constant 

Wavelength 

Laplace-transform operator 
electronic mass 

order of highest term kept in Fourier-Hermite expansion of f 
electron density; Fourier- Laplace transform of electron density 



n 0 mean electron density 

s complex argument of Laplace transform 

T absolute temperature 

t time 

v electron velocity 

v.. Q first component of normalized eigenvector corresponding to jth 

eigenvelocity Wj 

w phase velocity, w/k 

w n nth eigenvelocity 

x spatial coordinate 

a,j3 Fourier and Hermite indices, respectively 

-y imaginary part of w (y is negative for stability) 

AB arbitrary small increment in B 

AD n change in D n when w becomes w + Aw 

Aw arbitrary small change in w 

6 small positive quantity 

6 o!j3 Kronecker delta 

e plasma dielectric constant 

e(k,w=0); permittivity of vacuum 

X eigenvelocity 

p complex argument, w/k/2 
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<P = n(k,s) 

i// Fourier transform of time variation of 

\p. . . ordinate of histogram between w- and 

co argument of D n ; complex eigenvalue (root, pole, eigenfrequency); s/i 

u> n nth order approximation to eigenvalue 

( } average over velocity space 

Subscripts: 

i imaginary part 

r real part 

Superscript: 

* complex conjugate 


n 


w. 


j+1 


ANALYSIS 

Fourier-Hermite Representation 

Consider the ordinary one-dimensional Boltzmann equation (in natural units) for the 
single -particle distribution function f oftheelectrons: 

i +V i- E lk= C ® (1) 

where the right-hand side denotes a collision term that is a functional of f. The natural 
units of equation (1) are 


€ -0^L 

'n n e 2 


Debye length 



thermal speed 
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inverse plasma frequency 


g£L 

V n o e2 

The Debye length is a measure of the distance at which an introduced charge is screened, 
and the plasma frequency measures the rapidity of oscillations in electron density. A 
neutralizing uniform positive background is assumed. The external field E acting on 
the electrons is the self-field of the electrons themselves given by Poisson's equation: 


9E 

9x 



1 


( 2 ) 


If the number of particles in a Debye cube is large, as in a thermonuclear plasma, the 
collision term C may be neglected. In the case of identically vanishing C, equa- 
tions (1) and (2) are known as Vlasov equations. 

The Fourier-Hermite expansion of f is 


OO OC 


\j2iri = ^ ^ a^ exp(id!kx)H i g(v)exp(-v2/2) 

a=-°° /3=0 


where 


Hp(v) «= (-l)^exp(v2/2)^j exp(-v 2 /2) 
is Rodrigues' formula and 


(3) 


pCO 

- OO 


Hq, (v) H |3 ( v) exp (- v2/2) dv = n'. 


Substitution of equation (3) into equations (1) and (2) yields, after multiplication by 
exp(-io!kx)H^(v) and integration over x-v space, an infinite set of ordinary differen- 
tial equations for the coefficients of the Fourier-Hermite expansion. This set of 

equations for the matrix elements a^ is the Fourier-Hermite representation of the 
Vlasov equations (1) and (2). Truncation after the second row yields the Fourier-Hermite 
linearization. Further truncation at /3 = N yields the finite approximation to the linear- 
ized case. Finally, assumption of an exp(io)t) variation of the elements a^ leads to 
a dispersion relation of order (N + 1) . 

With Maxwell's distribution as background, the Fourier-Hermite dispersion rela- 
tion (from ref. 5 with a sign change on B) for cutoff of the expansion at order N is 
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det 


O) 

( l + k 2) 1/2 
0 
0 
0 


(l + k2) 1/2 

co + iB 

k \]2 

0 

0 


0 

k\/2 

co + 2iB 


0 

0 


kN 1 / 2 


0 

0 

kN 1 / 2 
co + NiB 


= 0 


( 4 ) 


In equation (4), co is an eigenfrequency (eigenvalue) and k is the wave number of an 
assumed small cos(kx) disturbance in electron density. If co is a root of equation (4), 
then so is -co*. The parameter B = 0 is a measure of the strength of the Fokker- 
Planck collision term (ref. 6), 


f, = B 


iL (vf ) + < v2> & 


lev'"' ' ' 9 v 2J 

which modifies the one-dimensional Vlasov equations (1) and (2) into 


9f df tt 8f f 
— + v-r £ --Et i - = f t 

0|; Sv I 


9x 


9v 


aE = r 

9x J_ 


f dv 


1 


The Vlasov case occurs at B = 0. When B = 0 the eigenvalues of equation (4) are all 
real, which suggests the equivalence of the Fourier-Hermite representation (ref. 5) to 
that of Van Kampen (ref. 2) rather than to that of Landau (ref. 1). However, as B 
increases from zero, the surprising behavior first reported in reference 3 and shown in 
figure 1 occurs. All but one of the poles (eigenvalues) move rapidly into the yco r -plane 
where a> = o; r - iy and f cc exp(icot). The anomalous pole is detained in the region of 
the Landau pole. Detailed investigation of the behavior of the anomalous pole for Fourier- 
Hermite expansions up to order 63 suggested that the equivalence of Van Kampen and 
Landau treatments was being observed within the Fourier-Hermite representation. A 
more powerful numerical technique has verified the behavior to orders over 1000, which 
is sufficiently near infinity to establish the nonuniformity convergent behavior of the 
approximation curve as N — oo and B 0. 


Numerical Algorithm 

The improved technique is based on recursion (ref. 7). H the determinant of order 
(N + 1) in equation (4) is denoted as D^ + -^, expansion on the last column yields 
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D N+1 = (W + iNB)D N - Nk 2 D N _! 


(5) 


In particular, 

Dj = oj 

D 2 = a)(w + iB) - (i + k 2 ) 


In principle, the value of D n+ j for arbitrary large N is obtained in about N steps, 
each of simple multiplication and addition. In contrast, the full matrix treatment of ref- 
erence 3 requires N 2 words just to record the elements of the matrix. 

To apply recursion equation (5) to the eigenvalue problem, the Fourier-Hermite 
root (eigenvalue) nearest in value to the real part of the Landau pole is identified for 
B = 0. A small positive increment in B, AB, and a small increment in a>, Au>, are 
introduced and the approximate derivative AD n+1/ Aw is computed. Use of the approx- 
imate derivative in the Newton- Raphson method gives a first approximation to the shifted 
eigenvalue. In general, 


-n + l = -n- D N + l( B+ ^-n)^; 

> 

AD N+1 = D N+l( B+AB ’ W n +Aw ) - D N+l( B+AB ’ w n) 


( 6 ) 


After the first step, iteration on the quasi- Newton-Raphson equations (6) yields a sharp 
value of u>(AB). A second increment of AB leads on to w(2AB), and in like manner 
the progress of gj(B) = u> T - iy on the yu> r -plane is tracked. The chief difficulties are 
found in scaling the large numbers which arise for large N and in the necessity of 
iteration to sharpen the values of gj(B). A sample FORTRAN program is shown in 
appendix A. 

Use of the recursion equation (5) simplifies the task of finding eigenvalues for the 
B = 0 (collisionless) case of equation (4). Easy evaluation of D n+ i(oj) at any N per- 
mits a movement from o> = 0 in steps of Aco > 0. Sign changes of D n+ j(co) localize 
each root within Aco, and iteration sharpens the root. The tridiagonal form of the matrix 
in equation (4) makes trivial (for k = 0.5) the calculation of the corresponding eigenvec- 
tors. In reference 5, histograms corresponding to an initial cos(kx) spatial disturbance 
of plasma density were constructed from the eigenvectors of the matrix of equation (4). 
Full matrix methods were barely adequate to yield histograms for the case of truncation 
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II 


at N = 99, and then only because B = 0 produces real eigenvalues and eigenvectors. 
With the recursion methods, histograms for N > 10 2 have now been constructed. 

The oscillations in density which follow an initial cos(kx) spatial density pertur- 
bation may be represented exactly by the following Fourier integral over a real phase 
velocity w: 



i//(w)exp(ikwt)dw 


(7) 


where kw = co, a real frequency. The reality of the Fourier-Hermite eigenfrequencies 
for B = 0 permits the construction of a histogram analogous to \p( w) in equation (7). 
In the limit of infinite N the histogram approaches w) . In appendix B the derivation 
of \p is given. As shown in reference 5, the ordinate of the histogram between the jth 
and (j + l)th eigenvelocities may be assigned as 


*4 


_ |( v j 2 ,o + v j>i,o) 


j,j+l 


,0 + v j+l> 

w. 1 - w. 
J+l J 


( w j+l > w j > 3 = 0,1,2, . . ., Nj 


where v n q is the first component of the normalized eigenvector corresponding to the 
nth eigenvelocity w n found by solving 


det 


(l + k-2) 
0 
0 
0 


w 

.,yi/2 


(l + k-2) 


- 2^/2 


w 

V/2 

0 

0 


0 

w 


0 

0 


N 1 / 2 


0 

0 

0 

N i/2 

w 


= 0 


( 8 ) 


Equation (8) is derived by similarity transformations of the matrix of equation (4) and by 
the introduction of the phase velocity w = w/k. 


RESULTS AND DISCUSSION 


The behavior of the Fourier-Hermite poles with vanishing collisions is indicated in 
figures 2 and 3 for approximations up to N = 1023. In figure 2 the behavior of the real 
part of the eigenfrequency is shown. A scale to accommodate all the values of N has 
been used in figure 2(a). The extreme deviations by which the curves reach Fourier- 


V 
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Hermite values at N = 255 and 511 as B goes to zero are quite striking. The curve 
for N = 1023 lies so nearly on a straight line to the Landau value that it would not show 
in figure 2(a) . Therefore, the behavior for the three highest values of N is shown at 
larger scale in figure 2(b). For N = 1023 it can be seen that the chance closeness of 
the Fourier-Hermite root to the Landau value causes a nearly straight approach to the 
Landau value. Actually, a deviation from a straight approach starts at about B = 0.0006. 

In figure 3 the behavior of the damping is shown for vanishing collisions. Again 
the two parts of the figure are small-scale and large-scale graphs. The curves for the 
higher N can be clearly seen only in the large-scale graph of figure 3(b). In every 
case the sharp swerve toward y = 0 from a straight approach to the Landau value is 
striking. 

It is clear from figures 2 and 3 that, numerically, the connection between the 
Van Kampen and Landau representations has been demonstrated within the Fourier- 
Hermite representation. Moreover, the pathological character of the curves u> r (B) 
and y(B) for large values of N as B goes to zero clearly indicates the nonuni- 
form ly convergent behavior of the approximation curve in the limit of infinite N at 
vanishing B. 

In figure 4 a histogram for N = 1023 is shown along with one for N s 99 that 
was presented in reference 5. The curve is symmetrical about X = 0, so only the posi- 
tive half is shown. Although a better representation of the curve for N = <» is obtained 
with N = 1023 than with N = 99, the peak is still rather crudely modeled despite a ten- 
fold increase in the number of terms in the Fourier-Hermite expansion. The root sepa- 
ration diminishes so slowly with increasing N that most roots can make only negligible 
contributions to the peak. 

The advantage of a collision term barely large enough to move a root very close to 
the Landau pole is suggested by figure 4. With a small collision term, a single Fourier- 
Hermite root can well represent the Landau damping, which is the most prominent feature 
of the behavior of the plasma at long times for weak initial disturbances of the plasma. 

CONCLUDING REMARK 

In addition to yielding conclusive indications of the nonuniformly convergent behavior 
of the Fourier-Hermite expansion as order increases to infinity and collisions decrease, 
the results have provided a guide to the choice of the size of collision term that best 
simulates the secular behavior of the plasma. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., March 23, 1972. 
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APPENDIX A 


CALCULATION OF EIGENVALUES 


The FORTRAN IV computer program used to compute the curves for N = 1023 in 
figures 2 and 3 follows. Although the dimensions 100 and 1024 have been used respec- 
tively for arrays APP (corresponding to successive root approximations) and DET (cor- 
responding to the successive determinants in the recursive process), these locations are 
a luxury for diagnostic purposes. Several locations would suffice if the successive root 
approximations and higher order determinants were written over the earlier values. For 
the case shown, central processor time on the CDC 6600 computer was less than 
30 seconds. 

PROGRAM CHASE ( I NPlJT * OUTPUT ) * 

♦ TRACKS F-H ROOTS INTO COMPLEX PLANE AS B INCREASES FROM ZERO * 

COMPLEX D * OMEGA * DERD * DEL TOM , DELMEG , APP * DET , OMEGA 0 * 

COMPLEX DFRN 1 * DERN2 * 

REAL K*K^OR, 

COMMON KSOR • B , ! S * SCALE * 

D I MFnS ION APP ( 1 00 ) t OFT (1024)* 

K = ^ • cr — O 1 , 

*V ALUF OF N AS USED BELOW CORRESPONDS TO ( N- 1 ) HFRM I TE FUNCTIONS* 

N = 1 0?4 * 

♦LANDAU POLE FOR K = 0.5 IS AT ( 1 • 4 1 566 * . 1 5336 > ♦ 

♦ F OUR I ER-Hf^RM I TE EIGENVALUE NEAREST REAL PART OF LANDAU POLE* 

OMEGA 0= <1.41 424 ,0* > * 

PRINT 900*M,!<* 

♦MAXIMUM NUMBER OF ITERATIONS ALLOWED TO FIX OMEGA, 

NOIT=100. 

♦ SHARPNESS CRITERION FOR OMEGA * 

CVLIm=2.F-oq, 

♦STEP SIZE IN APPROXIMATE nERIVATIVE* 

DFLT OM = < 0 • * 1 .E-06 ) * 

♦5— RANGE AND STFP SIZE* 

PM IN=* O002 * 

DFLTn = • 0OO02 * 

RMAXt: .0013, 

KSOR = K*K • 

OMFGA=OMFGAn, 
r = RmtN-DELTR * 

♦L OOP ON R* 

1 B=B+nELTR * 

♦NEWTON-RAPHSON LOOP* 

D02T = 1 * NO IT* 

*Uf2 * 

DERN2=D (N*OMEGA+DELTOM,DET ) * 

♦NUMRER OF SCALINGS* 

I S 2= I S * 

♦W 1 * 

DFRN f = D (N* OMEGA * OFT ) * 

TS1=IS* 
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APPENDIX A - Concluded 


*CHFC,K TO INSURE SAME SCALING ON W2,W1, 

IF( IS2.EQ. IS1 ) GO TO 4, 

IF< IS2-IS1 .GT.O > GO TO 6* 

PFRN?=DERN2*SCALF, 

GO To 4 * 

6 OFRN1 =derni*scale, 

^APPROXIMATE DERIVATIVE (W?-W1 )/DELTOM, 

4 DFRD= (DERN2-DFRN1 )/DeLTCM, 

3-NFWTON-RAPHSON ROOT I NCRFmFNT • 

DFLMeG=DFRN 1 /DEPO , 

OMFG/v ~ OMEGA - DEL MEG * 

APP ( 1 )=OMFGA , 

♦CHECK ON IMPROVEMENT OF OmEGA, 

2 IF (CABS (DELMEG ) .LT.CVL IM ) GO TO 3# 

G PRINT 1 00 , F3 , I , OMEGA , 

PRINT lOliISf 

100 FORMAT ( 3H B = F 1 0 . 6 , EX 1 1 H I T ERAT I ONS = I 3 * 5X6H0MEGA =F 1 1 • 8 , F 1 1 • 8 ) , 

101 FORMAT ( 1 X3H I S = I 3 ) , 

9CP FORMAT (//3H N= I 4/ 3H ^ = F6.3//)i 
IF (R,LT • BMAX — DEL TB/2 , ) , 

?GO TO 1,, 

STOP, 

FKjn , 

♦SUBPROGRAM FOR DISPERSION DETERMINANT, 

COMPLEX ^UNCTION D ( I ORD , OMEG A , DET ) , 

COMPLEX OMEGA , DET • 

RFAl ksor, 

COMMON KSQR,B, IS, SCALE, 

DIMENSION OFT (I ORD), 

*F IRST DETERMINANTS* 

DET ( 1 ) = OMEGA , 

DET ( G ) = OMEGA* (OMEGA+CMPLX (0# * ! • 0 ) *B ) - ( 1 *+KSQR ) * 

SCALF=1 • F — 18, 

IS=0, 

*RFCURSION LOOP, 

DO 1 1=3, I ORD , 

DET ( I ) =DET ( 1-1 )* ( OMEGA -f-FLO AT ( 1-1 >*CMPLX ( 0 • , 1 • ) *B ) -FLOAT ( 1-1 )* , 
2KSQR*DFT ( 1-2 ) , 

*CHFCK ON eIZF OF D, 

IF (CABS ( DET ( I ) )*LT,1 # E18) GO TO 1, 

♦SCALING Oe D, 

2 D^T ( T ) = DFT ( ! )*SCALF, 

DET ( T-l ) =DET ( 1-1 )*SCALE, 

♦SCALING COUNT, 

I S = ! s + 1 , 

1 CONTINUE, 

D=det( I ORD) , 

RETURN, 

END, 
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APPENDIX B 


DERIVATION OF xp(X) 

The linearized one-dimensional Vlasov equations, in natural units, for the electronic 
perturbation f on a uniform ionic background distribution are 


M + v |£-- ef' =o 
at 9x 

(Bla) 

/-»oo 

-|E = r f dx 

(Bib) 

J_oo 

where F is the unperturbed velocity distribution of the electrons, 
transform of f(x,v,t) may be written as 

The Fourier -Laplace 

f(k,s) = f dt f dx e _st e +lkx f(x,v,t) 

277 JO J-oo 



where the same symbol f is used for both function and transform. The corresponding 
inverse is 

P+ioo+s P<» 

f(x,v,t) = \ ds \ dk e +st e _lkx f(k,s) 

^-ioo+5 J-oo 

where the integral on s is known as Bromwich's integral. On application of the Fourier- 
Laplace transform to equations (Bl), the result is 

s f(k,s) - ikv f(k,s) - F’ E(k,s) = f Q (k,v) (B2a) 


pOO 

+ik E(k,s) = \ f(k,s)dv * n(k,s) 

^ — an 


To obtain equations (B2) the usual relations 


t 



s/{f} - f(x,v,0) 


9 



= -iky{f} 


(B2b) 


I 
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APPENDIX B - Continued 


are employed. The notation f Q (k,v) refers to ;?-(f(x,v,0))> , the Fourier transform of the 
initial conditions. In equations (B2), f(k,s) has an implicit dependence on velocity v, 
but n(k,s) and E(k,s) do not, since only in equation (B2b) has velocity been integrated 
upon. 

If one first solves for f(k,s) in equation (B2a), then integrates over v to find 
n(k,s) from equation (B2b), and, finally, substitutes s = ico, the result is 


n(k,s) 


1 f 0 (k,v)dv 
I(k,s) _ i J_ TO co - kv 

e(k,s) x + 1 r°° F^v 
k J_oo co - kv 


(s = ico) 


where the denominator e(k,s) is the so-called plasma dielectric constant. Under 
an initial cos(kx) density perturbation the assumed Maxwellian electronic velocity 
distribution 


F = -L= exp(-v2/2) 

|/2ir 

is untouched. Thus 

£ 0 (k,v) = F = - El 
and in terms of F' , 


n(k,s) 



1 -F 
v co - kv 



F f dv 
co - kv 


In the integral of the numerator of n(k,s) the coefficient of -F T may be written as 
v co - kv - kv co - kv/ 

so that n(k,s) may be written in terms of e(k,co) and e(k,0): 

nfk s) - ik 2 e ( k > a? ) ~ £ ( k >°) _ ik 2 6 ~ e 0 
’ co e(k,co) co e 
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APPENDIX B — Continued 


By means of a change in variable w = v/\j2, and the substitution 
w _ p _ 1 

jU ~ W fl - W 

the dielectric constant e(k,u>) may be expressed in terms of the complex function 

2 

Z(ju) = ~ f \ ~ /lm(p) <0; fj. = = x + iy\ 

pj.oo t - II I k/2 J 

as 

k 2 e(k ,p) = 1 + k 2 + iiZ(ii) 

In terms of Z = Z r + iZ^, 

k 2 e r = 1 + k 2 + xZ r - yZ^ 

k 2 e i = yZ r + xZ i 

The function Z is related to the complementary error function by 
Z(z) = ^e -z ^erfc(-iz) 

On the real axis of the p -plane, y = 0; so for real A = u>/k, 


k 2 e r (k,a;) = 1 + k 2 + -y= Z r /-pr\ 

V 2 IV 2 / 

(B3a) 

k2£ i (k ’“ , = # z #) 

(B3b) 

k 2 e r (k,0) = 1 + k 2 

(B3c) 

k 2 Cj(k,0) = 0 

(B3d) 


Inversion by Bromwich's integral with the contour along the imaginary axis of the s-plane 
(5 = 0) yields 
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APPENDIX B - Concluded 


n(k,t) =^r e i W t dw = s£ f 

2» J.„ “* 2i J_, 




U> 


(B4) 


At all t < 0, n(k,t) = 0; thus 


poo (p + up 

0 = l — e _la,t du) 

J_oo W 


(t>0) 


Taking the complex conjugate yields 


r 

J-on W 


(t >0) 


Subtraction of ik^/27T times this null member from equation (B4) for n(k,t) yields 


2 r°° <p • ( co) . . n°° 

n(k,t) = - — \ — i- — e la, ^da> = \ ^/(w)e lut dw 

^ - OO — 0O 


Writing in terms of the values of e r and e i given in equation (B3) yields 


i//(o>) = 


1 + k 


- 2 Zi ikh) i 


i r\j2 e^(k,u>) + e2(k,cu) k 


where the function Z + = -Z* is one available on the computer library tape at Langley 
Research Center. The substitution w = kX (using the same symbol for \p) yields 


n(k,t) = f 4/(x)e ik ^ t d\ 

•J -ca 


where 


M a ) = 


_ 1 H- k~ 2 Z i(^) 


Tl\j 2 


r i 


is the function corresponding to infinite N that is plotted in figure 4 . 
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Figure 1.- Movement of Fourier-Hermite poles 
into complex frequency plane as collisions 
increase. N = 63; k = 0.5. 










(a) N = 99. 



A 

(b) N = 1023. 

Figure 4.- Exact Fourier transform compared with 
corresponding Fourier-Hermite histograms for 
initial cos(kx) density disturbance, k = 0.5. 
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